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LONG-TERM  GOALS 

Our  long-tenn  objective  is  to  establish  a  database  of  sediment  and  bottom  properties  including 
sediment  compressional  wave  speed,  shear  wave  speed,  attenuation  and  density.  Using  this  we  propose 
to  provide  the  best  sediment  model  for  any  region  with  estimates  of  model  uncertainty.  These  data 
would  be  appropriate  for  incorporation  in  state-of-the-art  propagation  models  for  acoustic  effects 
prediction  on  the  marine  environment. 

OBJECTIVES 

1 .  Develop  a  surficial  sediment  model  for  the  world’s  oceans  based  on  optimal  estimation 
techniques,  geological  models  and  data. 

2.  Integrate  the  sediment  model  into  the  ESME  module. 

3.  Provide  estimates  for  model  uncertainty. 

4.  Apply  this  model  in  selected  regions  to  test  its  accuracy  and  effectiveness. 

APPROACH 

As  a  first  step  we  collect  all  the  available  sediment  property  data  at  the  region  of  interest.  This  includes 
all  the  geoacoustic  (compressional  wave  speed,  shear  wave  speed,  attenuation  etc.)  and  geotechnical 
properties  (bulk  density,  permeability,  porosity  etc.).  In  addition  other  useful  sources  of  information 
like  sediment  thickness  maps,  geophysical  survey  records  etc.  are  also  used.  The  geotechnical 
properties  can  be  used  to  derive  geoacoustic  properties  using  a  suitable  model.  These  data  can  be 
obtained  from  literature  or  from  available  sediment  property  databases.  All  this  data  is  then  used  to 
efficiently  estimate  the  sediment  properties  at  any  location  within  the  region  of  interest. 

WORK  COMPLETED 

Sediment  models  suitable  for  the  Southern  California  Bight  (SOCAL)  region  and  the  Middle  Atlantic 
Bight  region  were  developed  and  integrated  into  the  ESME  module  as  part  of  the  two  test  cases.  The 
third  test  case  was  at  the  VAST  site.  A  program  was  written  to  generate  sediment  parameters  for  the 
VAST  test  region.  The  location  lies  between  latitudes  33  N  and  40  N  and  71  W  and  77  W.  The 
program  is  written  in  such  a  way  that  given  the  location  and  the  frequency  of  interest,  the  sediment 
properties  are  estimated  using  all  the  available  data. 

In  addition,  an  EOF  based  representation  was  developed  for  the  East  China  Sea  region.  Data  for  this 
came  from  the  ASIAEX  experiment  sponsored  by  ONR.  Sediment  properties  were  interpolated  using 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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this  approach  and  compared  with  results  from  long  range  sediment  tomography.  The  results  are 
described  in  detail  in  the  following  sections. 


RESULTS 

A.  Sediment  models  for  ESME  test  region  (VAST) 

Sediment  models  suitable  for  the  Southern  California  Bight  (SOCAL)  region  and  the  Middle  Atlantic 
Bight  region  were  developed  and  integrated  into  the  ESME  module  as  part  of  the  two  test  cases.  The 
third  test  case  was  at  the  VAST  site.  A  program  was  written  to  generate  sediment  parameters  for  the 
VAST  test  region.  The  location  lies  between  latitudes  33  N  and  40  N  and  71  W  and  77  W.  The 
program  is  written  in  such  a  way  that  given  the  location  and  the  frequency  of  interest,  the  sediment 
properties  are  estimated  using  all  the  available  data. 

The  bottom  properties  are  specified  for  a  simple  two  layer  model.  Some  of  the  properties 
(compressional  and  shear  speeds  in  the  basement)  are  assumed  constant  throughout  the  area.  The 
sediment  information  available  include  core  data,  sediment  maps,  sediment  type  infonnation  etc.  For 
this  location  we  were  able  to  find  some  historical  data  which  are  listed  below. 

1 .  Core  data: 

Atlantic  Margin  Coring  Sites  6006  to  6012  and  6020  to  6021  ’  .  The  sediment  type  information, 
sediment  properties,  layer  thickness  etc.  were  obtained  from  these. 

Deep  Sea  Drilling  sites,  902  to  906. 

2.  Sediment  type  information:  from  USGS  Open  File  Report3. 

3.  Sediment  sound  velocity  functions  for  North  Atlantic  (Houtz)4 

4.  Thickness  of  surficial  sand  sheet  (Knebel  and  Spiker)5 

5.  Sediment  and  geoacoustic  models: 

Hamilton  model6:  Used  for  predicting  the  compressional  and  shear  wave  speeds  as  well  as 
attenuation  as  function  of  frequency  for  different  sediment  types. 

Bowles7:  Used  for  predicting  attenuation  as  a  function  of  frequency  and  shear  speeds  for 
known  sediment  types  . 

Biot-  Stoll  model8  to  compute  the  compressional  and  shear  wave  speeds  from  sediment 
paramaters. 

6.  Sediment  layer  thickness  maps:  Used  for  predicting  the  sediment  thickness. 


Sediment 


Compr.  Speed,  Shear  Speed 
Attenuation,  Density 


Basement 


Compr.  Speed,  Shear  Speed 
Attenuation,  Density 


Figure  1.  Two  layer  model  to  represent  the  sediment  properties  used  in  this  study. 

The  program  uses  these  data  to  arrive  at  the  sediment  parameters  at  a  given  location.  The  sediment 
compressional  speed  is  estimated  using  the  core  data,  sediment  type  infonnation  and  the  Hamilton’s 
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values  for  these  sediment  types.  Sediment  thickness  is  estimated  based  on  the  core  data  and  sediment 
thickness  map.  Shear  speed,  and  Compressional  and  shear  attenuation  are  determined  mainly  using  the 
sediment  type  information  and  predictions  by  Hamilton  and  Stoll.  The  bottom  properties  are  specified 
for  a  simple  two  layer  model.  Some  of  the  properties  (compressional  and  shear  speeds  in  the  basement) 
are  assumed  constant  throughout  the  area.  The  output  parameters  are  written  in  a  text  file. 
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Figure  2.  Left  panel  shows  the  location  of  the  ESME  test  case  and  the  bathymetry.  Right  panel 
shows  the  sediment  sound  speed  obtained  using  the  sediment  model  developed. 

B.  Empirical  Orthogonal  Function  (EOF)  based  representation  of  sediment  properties. 

EOFs  have  been  used  to  represent  variations  in  ocean  sound  speeds  for  a  long  time.  They  have  also 
been  used  to  represent  variations  in  sediment  properties  recently9.  EOFs  provide  an  efficient  means  of 
representing  measurements  as  a  set  of  components  or  modes  that  are  uncorrelated.  Observations  of  any 
field  can  be  expressed  as  a  sum  of  the  mean  and  a  perturbation.  These  field  perturbations  can  then  be 
expressed  as  a  sum  over  the  EOF  basis  vectors  incorporating  sufficient  number  of  modes  to  capture 
most  of  the  variations.  The  procedure  for  implementing  this  in  the  context  of  sediment  modeling  is 
described  in  detail  by  Jay  et  al.9  We  have  tried  to  use  this  approach  to  East  China  Sea  location  which 
was  the  site  for  the  ASIAEX  experiment  sponsored  by  ONR.  As  part  of  the  experiment  various 
geoacoustic  measurements  were  taken  which  included  gravity  and  piston  cores,  water  gun  and  chirp 
surveys. 

Figure  3  shows  the  result  of  the  EOF  based  interpolation  and  its  comparison  with  sediment 
tomography  inversions.  The  core  profiles  from  locations  indicated  in  the  left  panel  of  the  figure  are 
used  to  obtain  the  EOFs.  The  sediment  speeds  at  three  locations  along  a  source-  receiver  path  is  then 
estimated  using  the  EOFs  using  the  approach  given  in  Jay  et  al9.  The  right  panel  shows  the  comparison 
of  these  interpolated  profiles  with  the  tomography  inversions.  The  tomography  inversions  produce  a 
mean  profile  calculated  from  the  Genetic  Algorithm  solutions  and  also  a  profile  from  the  hybrid 
approach.  Potty  et.al10.  have  discussed  this  inversion  scheme  in  detail. 
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•  ECS  Shot  #  41 


compressional  wave  speed  (m/s) 


Figure  3.  Left  panels  shows  the  distribution  of  the  cores  in  the  experimental  area.  The  radius  of  the 
circle  is  approximately  30  km.  The  blue  asterisks  are  explosive  charges  used  for  sediment 
tomography.  The  EOF  based  approach  is  used  to  interpolate  the  sediment  sound  speed  at  three 
locations  (black  diamonds)  along  the  propagation  path  from  one  of  the  explosive  charge  to  the 
receiver  which  was  at  the  center  of  the  circle.  Right  panel  shows  the  interpolated  sediment  speed 

data  compared  with  tomography  inversions 


C.  Future  work  planned 

Our  future  goals  include  the  following: 

1 .  Use  the  EOF  based  representations  to  efficiently  incorporate  geophysical  data  into  the  sediment 
model 

2.  Provide  some  measure  of  uncertainty  to  the  estimates. 

IMPACT/APPLICATIONS 

Our  effort  is  intended  to  provide  the  best  possible  estimates  of  sediment  data  needed  for  the 
propagation  modeling  component  of  the  ESME  initiative.  This  will  enable  the  propagation  modelers  to 
use  ‘best  guesses’  when  direct  estimates  of  the  sediment  data  is  not  readily  available. 

TRANSITIONS 

We  expect  that  this  sediment  model,  when  complete,  will  be  useful  to  the  acoustic  community  as  a 
whole  as  an  important  database  of  sediment  property  information. 

RELATED  PROJECTS 

We  have  used  data  from  the  ASIAEX  experiment  sponsored  by  ONR  to  test  the  EOF  based 
representations  for  sediment  properties.  These  representations  will  eventually  be  useful  for  geoacoustic 
modeling  of  sediments. 
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